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Abstract 

The  small  polymer  electrolyte  fuel  cell  (PEFC)  was  experimented  by  changing  operation  temperature  and  Ho  and  Oi  concentration  in  supply 
gas,  and  the  PEFC  reaction  characteristic  model  that  can  express  these  influences  was  created.  Moreover,  PEFC  reaction  and  flow  analysis 
model  was  made  with  this  reaction  characteristic  model  and  thermal  flow  analysis  model.  Furthermore,  in  order  to  improve  the  performance 
and  the  safety  of  PEFC,  five  kinds  of  separators  were  evaluated  from  the  viewpoint  of  state  of  gas  flow,  uniformity  of  current  density  and 
temperature,  reduction  of  pressure  loss  and  exhaustibility  of  liquid  water  with  this  PEFC  reaction  and  flow  analysis  model.  In  case  of  parallel 
separator,  flow  rate  and  current  density  distribution  was  not  uniform,  and  exhaustibility  of  the  liquid  water  was  low.  In  case  of  serpentine 
separator,  pressure  drop  was  large.  In  addition,  semi-serpentine  separator  was  designed  for  the  shape  which  combines  parallel  separator  and 
serpentine  separator,  and  the  temperature  and  current  density  distribution  of  this  separator  was  more  uniform  than  any  other  separators,  so 
this  separator  was  excellent. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

At  present,  humankind  have  serious  problems  of  environ¬ 
ment,  such  as  global  warming  and  acid  rain,  and  of  lack  of 
fuel  sources,  such  as  petroleum  and  natural  gases.  In  order 
to  contribute  to  the  solution  of  these  problems,  fuel  cell  is 
expected  to  be  practical  use  because  it  has  low  emission 
of  environmental  pollutant  and  high  conversion  efficiency 
from  chemical  energy  to  electrical  energy.  Especially,  poly¬ 
mer  electrolyte  fuel  cell  (PEFC)  is  developed  as  power  source 
of  transportation  and  stationary  family  power  supply  in  vari¬ 
ous  fields,  because  it  has  low  operation  temperature  and  high 
power  density.  The  materials  development  and  the  configu¬ 
ration  design  of  separator,  which  is  component  of  PEFC,  are 
also  researched. 

In  order  to  improve  the  power  generation  performance, 
the  durability  and  the  safety  of  PEFC,  it  is  important  to 
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design  the  best  separator  shape  to  satisfy  the  following  three 
conditions  concerning  the  whole  PEFC  system:  (1)  the  cur¬ 
rent  density  distribution  and  the  gas  flow  rate  distribution 
can  be  unified;  (2)  liquid  water  can  be  removed  smoothly; 
(3)  pressure  drop  of  supply  gas  can  be  reduced.  However, 
there  are  very  few  researches  that  examine  the  design  of  gas 
flow  channel  by  considering  internal  phenomena  while  the 
electric  power  is  generated.  As  the  flow  and  mass  transfer 
phenomena  occur  complicatedly  between  microscopic  area 
and  macroscopic  area,  it  is  very  difficult  to  visualize  and 
measure  the  phenomenon.  Consequently,  numerical  analy¬ 
sis  is  useful  to  examine  it.  Bernardi  and  Verbrugge  [1,2] 
and  Springer  et  al.  [3]  developed  one-dimensional  model  to 
the  direction  of  membrane  thickness,  and  examined  concen¬ 
tration  distribution  and  water  management  in  PEFC.  Fuller 
and  Newman  [4]  analyzed  and  developed  two-dimensional 
model  to  the  direction  of  membrane  thickness  and  gas  flow 
channel.  Nguyen  and  White  [5],  and  Yi  and  Nguyen  [6]  devel¬ 
oped  heat  and  water  transport  models  (2-D)  that  accounted 
for  various  operation  conditions  and  membrane  hydration 
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Mj 

molecular  weight  of  species  j  (kg  mol-1) 

membrane  dry  equivalent  weight 

n 

number  of  electrons  participating  in  a  reaction 

11(1 

electro-osmotic  drag  coefficient 

p 

pressure  in  Eq.  (20)  (Pa) 

a  P 

water  ejection  differential  pressure  in  Eq.  (35) 
(Pa) 

p 

pressure  (Pa) 

7^0,  sat 

saturated  vapor  pressure  in  stream  (Pa) 

<71 

heat  flux  from  MEA  to  gas  phase  (J  m-2  s-1) 

<72 

heat  flux  from  back  plate  to  gas  phase 
(J  m-2  s-1) 

93 

heat  value  generated  by  reaction  per  unit  area 
(J  m-2  s-1) 

q4 

heat  flux  from  gas  phase  to  solid  phase 
(J  m-2  s-1) 

<75 

heat  flux  from  back  plate  to  solid  phase 
(J  m-2  s-1) 

<76 

latent  heat  value  of  condensation  per  unit  area 
(J  m-2  s-1) 

n 

molar  flux  of  species  j  (mol  m-2  s-1) 

R 

gas  constant  (8.314  J  mol-1  K-1) 

t 

time  (s) 

membrane  thickness  (m) 

T 

temperature  (K) 

]b 

cooling  plate  temperature  (K) 

T s 

gas  phase  temperature  (K) 

T 

solid  phase  temperature  (K) 

V 

flow  velocity  vector  (ms-1) 

V 

flow  velocity  (ms-1) 

V 

operation  voltage  (V) 

X 

horizontal  distance  (m) 

A.r 

distance  on  calculation  mesh  to  neighboring 
gas  channel  (m) 

y 

horizontal  distance  (m) 

z 

vertical  distance  (m) 

Greek  letters 

a 

net  water  transfer  coefficient 

aa 

transfer  coefficient  in  Eq.  (5) 

a\ 

transfer  coefficient  in  Eq.  (6) 

a2 

transfer  coefficient  in  Eq.  (11) 

*7  act 

activation  overvoltage  (V) 

Ocon 

concentration  overvoltage  (V) 

Oohm 

ohmic  resistance  overvoltage  (V) 

X 

water  content  in  the  membrane 

p 

viscosity  of  mixture  (Pa  s) 

p 

density  of  mixture  (kg  m-3) 

nm 

Pdiy 

membrane  dry  density  (kgm-3) 

crm 

membrane  conductivity  (1  <T2—  1  m-1) 

effective  membrane  conductivity  ( 1  !T2- 1  m- 1 ) 

r 

viscous  stress  (Pa) 

r 

viscous  stress  tensor  (Pa) 

§ 

activity  water  in  stream 

Nomenclature 

Ae 

effective  surface  area  per  unit  projection  area 
of  a  electrode 

effective  surface  area  per  unit  amount  of  plat¬ 
inum  (m2  g-1) 

bc 

condensation  rate  constant  (s—  1 ) 

Ca 

molar  concentration  of  species  A  (mol  m-3) 

C6 

molar  concentration  at  interface  between  GDL 
and  electrode  (molm-3) 

co2 

molar  concentration  of  oxygen  in  electrolyte 
(mol  m-3) 

^f.ref 

co2 

reference  molar  concentration  of  oxygen  in 
electrolyte  (molm-3) 

Cg 

molar  concentration  in  the  gas  phase 
(mol  m-3) 

rs 

co2 

molar  concentration  of  oxygen  in  the  gas  phase 
(mol  m-3) 

cp 

specific  heat  at  constant  pressure  (J  kg-1  K-1) 

r(m) 

L'H20 

molar  concentration  at  anode  or  cathode  side 
interface  of  membrane  (molm-3) 

c« 

parameter  used  in  the  expression  for  cathodic 
transfer  coefficient 

Ca 

fitting  parameter  in  Eq.  (7) 

D 

diffusion  coefficient  (m2  s-1) 

Da 

diffusion  coefficient  of  species  A  (m2  s-1) 

Dh 

hydraulic  diameter  (m) 

nm 

uh2o 

effective  diffusion  of  water  in  membrane 
(m2  s-1) 

E 

electromotive  force  (V) 

standard  electromotive  force  (V) 

Eah 

the  value  of  reduction  change  of  water  enthalpy 
to  voltage  (V) 

f 

effective  porosity 

F 

Faraday’s  constant  (96,485  C  mol-1) 

h 

heat  transfer  coefficient  (J  m-2  s-1  K- 1 ) 

A7/h2o 

change  of  water  enthalpy  between  vapor  and 
liquid  (J  mol-1) 

i 

current  density  (A  cm-2) 

it 

exchange  current  density  (A  cm-2) 

•ref+ 

l0 

exchange  current  density  at  reference  O2  con¬ 
centration  (A  cm-2) 

<L 

limit  current  density  (A  cm-2) 

k 

thermal  conductivity  (J  m-1  s-1  K-1) 

bio 

proportionally  constant  of  exchange  current 
density 

ka 

fitting  parameter  in  Eq.  (7) 

ld 

gas  diffusion  layer  thickness  (m) 

7g 

gas  channel  depth  (m) 

Is 

thickness  of  solid  phase  (m) 

L 

inlet  width  (m) 

fftpt 

amount  of  platinum  per  unit  electrode  area 
(gm-2) 
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Superscripts 

a 

anode 

b 

back  plate 

c 

cathode 

g 

gas  phase 

m 

membrane 

s 

solid  phase 

sep 

separator 

Subscripts 

A 

species  A 

j 

species  j 

h2 

hydrogen 

h2o(1) 

liquid  water 

H20(v) 

vapor  water 

MEA 

membrane  electrode  assembly 

n2 

nitrogen 

o2 

oxygen 

conditions.  On  the  other  hand,  it  is  thought  that  analysis 
with  the  computational  fluid  dynamics  (CFD)  technique  is 
important  in  order  to  calculate  the  transport  phenomena  in 
detail,  and  such  study  is  increasing  recently.  Um  et  al.  [7] 
and  Wang  et  al.  [8]  have  developed  two-dimensional  model 
with  CFD,  which  included  two-phase  flow.  Dutta  et  al.  [9] 
made  a  three-dimensional  computational  model  based  on  the 
commercial  software  package.  Fluent.  Berning  et  al.  [10] 
presented  a  non-isothermal,  three-dimensional  models  and 
calculated  the  distribution  of  current  density  and  concentra¬ 
tion  in  the  straight  channel.  Mazumder  and  Cole  [11]  exam¬ 
ined  the  liquid  water  transport  with  three-dimensional  model. 
Li  et  al.  [12]  analyzed  in  small  cell  with  three-dimensional 
analysis.  Um  and  Wang  [13]  compared  the  performance  of 
straight  flow  channel  with  that  of  interdigiteted  flow  channel 
by  three-dimensional  analysis.  Kulikovsky  [14]  analyzed  ser¬ 
pentine  channel  under  the  condition  that  gas  inlet  flow  rate  in 
each  channel  was  constant.  These  PEFC  numerical  analysis 
models  contributed  to  the  optimization  of  component  design 
and  operation  condition,  and  to  the  examination  of  issues 
included  in  present  cell.  However,  these  studies  calculated 
internal  phenomena  in  short  straight  gas  channel  and  one 
serpentine  channel  or  in  small  cell,  and  there  are  very  few 
researches  that  evaluated  various  channel  shapes  of  actual 
size  cell  under  considering  realistic  calculation  time  and  cal¬ 
culation  resource. 

In  this  study,  the  small  PEFC  was  experimented  with 
changing  operation  temperature,  humidify  temperature  and 
Ht  and  O2  concentration  in  supply  gas.  The  influence  to 
the  i—V  characteristic  was  considered,  and  the  PEFC  reac¬ 
tion  model  that  can  express  these  influences  was  created. 
Next,  the  actual  size  PEFC  reaction  and  flow  analysis  model 
was  made  with  this  reaction  model  and  thermal  flow  analysis 
code.  In  addition,  the  water  ejection  condition  was  investi¬ 
gated  by  the  one  channel  gas  flow  experiment.  Finally,  five 


kinds  of  separators  were  evaluated  with  this  PEFC  reaction 
and  flow  analysis  model  from  the  viewpoint  of  gas  flow 
condition,  current  density  distribution,  cell  temperature  dis¬ 
tribution,  pressure  drop,  ejection  performance  of  liquid  water. 
And  the  best  separator  shape  that  can  be  safety  and  high  per¬ 
formance  was  examined. 


2.  PEFC  reaction  model 

2.1.  Structure  of  small  PEFC  and  experimen  t 

Fig.  1  shows  the  structure  of  small  PEFC  in  this  experi¬ 
ment.  Membrane  electrode  assembly  (MEA),  which  is  coated 
with  Pt  electrode  on  both  sides  of  electrolyte  membrane,  is  in 
the  center  of  PEFC.  MEA  is  sandwiched  between  two  carbon 
gas  diffusion  layers  (GDL)  that  improve  repellency  and  gas 
diffusibility,  and  furthermore,  this  set  is  sandwiched  between 
two  separators  that  have  a  straight  groove;  a  depth  of  gas  flow 
groove  is  1  mm,  a  width  is  2  mm  and  a  length  is  50  mm.  Anode 
gas  (hydrogen  as  fuel)  and  cathode  gas  (oxygen  as  oxidant) 
flow  in  the  each  separator  channel,  respectively.  Specification 
of  MEA  and  GDL  is  shown  in  Table  1 .  Small  cell  was  used 
because  physical  quantities,  which  are  gas  flow  rate,  concen- 


Membrane  electrode 
assembly 

Gas  diffusion  layer 


Separator 


Fig.  1 .  Schematic  diagram  of  small  single  cell. 


Table  1 

Specification  of  MEA  and  GDL  of  small  cell 


MEA 

Size  of  membrane  (mm2)  50  x  100 

Thickness  of  membrane  (p.m  )  30 

Size  of  catalyst  layer  (mm2)  10x50 

Amount  of  Pt  (g  m-2)  3.0 

GDL 

Size  (mm2)  10x50 

Thickness  (p-m)  300 
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Fig.  2.  Schematic  diagram  of  experimental  apparatus. 


tration  and  temperature,  have  to  be  as  uniform  as  possible  in 
channel.  Fig.  2  shows  schematic  of  experimental  apparatus. 
Anode  gas  was  mixture  gas  of  hydrogen  and  nitrogen,  and 
hydrogen  compositions  were  3,  10,  30  and  100%.  Cathode 
gas  was  mixture  gas  of  oxygen  and  nitrogen,  and  oxygen 
compositions  were  5,  10  and  21%.  As  electrolyte  membrane 
is  affected  by  moisture  condition,  humidification  of  supply 
gases  is  necessary  in  order  to  keep  the  water  content  and 
the  ionic  conductivity  of  membrane.  Accordingly,  external 
humidifier  was  set  at  the  front  of  gas  inlet  of  cell,  and  cell  and 
humidifier  were  set  in  thermostat  and  water  bath,  respectively. 
The  cell  voltage  and  current  were  measured  with  multime¬ 
ter  and  load.  The  operating  condition  is  shown  in  Table  2. 
In  order  to  decrease  the  concentration  distribution  due  to  cell 
reaction,  hydrogen  and  oxygen  utilizations  were  set  to  be  less 
than  5%. 

2.2.  Assumptions  of  PEFC  reaction  model 

Fig.  3  shows  PEFC  reaction  simulation  model.  Operation 
voltage  is  calculated  by  following  equation: 

V  =  E  -  'fact  -  'fact  -  'Icon  “  'Icon  “  'Cm  (!) 

where  V  is  operating  voltage,  E  the  electromotive  force,  pa ct 
the  activation  overvoltage,  rjcoo  the  concentration  overvolt- 


Table  2 


Operation  conditions  of  small  cell  experiments 


Cell  temperature  (K) 

298,313,333,363 

Humidify  temperature  (K) 

298,313,333,363 

O2  composition  (%) 

5,  10,  21 

H2  composition  (%) 

3,  10,  30,  100 

Pressure  (MPa) 

0.1 

Anode  gas  flow  rate  (m3  s_  1 ) 

3.33  x  KT6 

Cathode  gas  flow  rate  (m3  s-1) 

8.35  x  KT6 

Fig.  3.  PEFC  reaction  model. 


age,  r/0|lm  the  resistance  overvoltage  and  the  superscript  a,  c 
and  m  express  anode,  cathode  and  membrane,  respectively. 
These  overvoltage  are  caused  by  electrode  reaction  and  trans¬ 
port  phenomena  at  the  parts  shown  in  Fig.  3,  respectively. 
In  this  study,  the  following  parameters  in  these  overvoltage 
model  derived  in  the  next  section  were  used  as  fitting  param¬ 
eter  to  help  the  calculation  data  agree  with  experimental  data 
in  various  conditions, 

f\f,a\a\,Ca,CaK,kio 

where  /a  is  the  effective  porosity  of  anode  side  GDL,  /c 
the  effective  porosity  of  cathode  side  GDL,  «a  the  transfer 
coefficient  of  anode  concentration  overvoltage,  o?|  the  trans¬ 
fer  coefficient  of  cathode  concentration  overvoltage,  Ca  the 
coefficient  to  calculate  the  transfer  coefficient  of  activation 
overvoltage,  C„  and  kn  the  correction  coefficients  to  calcu¬ 
late  the  membrane  ionic  conductivity  and  kj0  is  the  correction 
coefficients  to  calculate  the  exchange  current  density.  Though 
there  were  some  assumptions  and  models,  which  were  not 
strict  theoretically,  the  following  models  were  developed  in 
order  only  to  use  these  value  as  fitting  parameter. 

These  overvoltage  equations  are  derived  on  the  basis  of 
the  following  assumption. 

1 .  Gas  concentration  at  the  interface  between  gas  flow  chan¬ 
nel  and  GDL  is  equal  inlet  gas  concentration. 

2.  The  reaction  area  reduction  caused  by  flooding  of  elec¬ 
trode  is  ignored,  and  the  diffusion  prevention  caused  by 
water  condensation  is  also  ignored. 

3.  The  heat  of  reaction  is  discharged  out  of  PEFC,  and  the 
temperature  of  every  PEFC  component  is  constant. 

4.  In  membrane,  ionic  conductivity,  electro-osmosis  coeffi¬ 
cient  and  water  effective  diffusion  coefficient  that  depend 
on  membrane  humidity  are  determined  by  water  activity 
of  anode  side. 

5.  Anode  activation  overvoltage  is  ignored  because  it  is  lower 
than  any  other  overvoltage. 

6.  The  gas  crossover  is  disregarded. 
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2.3.  Derivations  ofPEFC  reaction  model 


modified  the  equation  of  Springer  et  al.  is  used: 


The  electromotive  force  is  shown  by  the  following  Nernst 
equation: 


CTem  =  (0.514-Co-A.  -  0.326 ka)  x  exp 


(7) 


E  —  E° 


RT 

nF 


In 


Pus-(PoJ/2 


(2) 


where  £°  is  standard  electromotive  force,  R  the  gas  constant, 
T  the  temperature,  n  the  number  of  electrons  participating  in 
a  reaction,  F  the  Faraday’s  constant,  the  anode  hydro¬ 
gen  partial  pressure  and  PPh  is  the  cathode  oxygen  partial 
pressure. 

Hydrogen  that  flows  along  the  anode  channel  moves  to 
interface  of  MEA  and  GDL  through  the  GDL,  and  hydrogen 
concentration  in  the  GDL  decreases  toward  MEA.  Conse¬ 
quently,  assuming  hydrogen  concentration  gradient  is  linear 
and  constant  in  GDL  with  thickness  /d,  current  density  can 
be  calculated  with  ,  which  is  concentration  on  interface 
of  channel  and  GDL,  and  Cjj  ,  which  is  concentration  on 
interface  of  GDL  and  MEA: 

i  =  2  Ffa  D\  i,  Chz  ~  Ch2  (3) 


where  D\\,  is  diffusion  coefficient  and  /a  is  the  effective 
porosity  of  anode  GDL.  When  is  zero,  concentration 
gradient  and  current  density  becomes  maximum.  The  anode 
limiting  current  density  is  shown  by  the  following  equation. 
Similarly,  cathode  limiting  current  density  can  be  calculated: 


cl 

IUH2)  =  2Ff-DUl  =  2FfDHl 


^Pk 

/d  RT 


Cs  1  P8 

iu  o2)  =  4F/cDo2^  =  4FfDo2--^. 


(4) 


where  C„  and  kr  are  constant.  C„  and  kn  were  equal  to  1 
in  the  equation  of  Springer  et  al.  However,  as  the  membrane 
that  was  used  in  Springer’s  experiment  was  different  from 
one  of  this  experiment,  so  the  ionic  conductivity  of  this  study 
is  calculated  with  these  parameter  that  were  fitting  parameter 
satisfied  the  coincidence  with  experimental  and  calculation 
results.  Nguyen  and  White  [5]  applied  water  content  of  anode 
side  to  water  content  in  Eq.  (7).  Water  content  of  anode  side 
is  calculated  with  the  anode  water  activity. 

Pa  =  0.043  +  17.8fa  -  39.8(£a)2  +  36.0(§a)3  (£a  <  1) 

Pa  =  14+  1.4(£a-  1)  (£a  >  1) 

(8) 


Anode  water  activity  is  shown  by  the  following  equation: 

(9) 

rH20,sat 


Saturated  vapor  pressure  is  calculated  by  Antoine  equation. 
Resistance  overvoltage  is  shown  by  the  following  equation: 

tm  . 

^ohm  =  — ~i  (10) 

IT™ 


where  tm  is  thickness  of  membrane. 

Anode  activation  overvoltage  is  calculated  by  the  follow¬ 
ing  Tafel  equation: 


7 act  — 


RT 

afF 


In 


A  i+ 


(ID 


where  Ae  is  effective  surface  area  per  unit  projection  area  of 
electrode. 


Anode  concentration  overvoltage  can  be  calculated  by  the 
following  equation  with  the  limiting  current  density: 


7 


a  _ 

con 


RT 

oP2F 


In 


!L(H2) 


(5) 


Ae  —  nZptAg  (12) 

where  mpt  is  the  amount  of  platinum  per  unit  electrode  area 
and  As  is  effective  surface  area  per  unit  amount  of  platinum. 
As  is  obtained  in  the  effective  surface  area  data  of  Marr  and  Li 
[15].  /q  is  oxygen  exchange  current  density,  and  it  is  obtained 
by  the  following  equation  in  reference  [16]: 


Similarly,  cathode  concentration  overvoltage  can  be  calcu¬ 
lated  by  the  following  equation: 


*7 


C 

con 


RT  f  i  \ 
— -  In  1 - 

“l2P  V  lU02)J 


(6) 


•ref+ 


c, 


o2 


C, 


rf.ref 
02  / 


ref 

lo2 


o2 


exp(14.1  -  666/P) 


03) 

(14) 


where  «a  and  aj  are  transfer  coefficient  of  anode  and  cathode, 
respectively. 

The  ionic  conductivity  of  electrolyte  membrane  is  calcu¬ 
lated  by  the  Springer  et  al.  method  [3],  In  this  study,  in  order 
to  agree  with  experimental  and  calculation  data,  the  following 


Cf02  =  0.07  C82  (15) 

f  ref 

where  C(h  is  the  reference  oxygen  molar  concentration  in 
the  electrolyte,  C q,  the  oxygen  molar  concentration  in  the 
electrolyte  and  Cf)n  is  the  molar  concentration  of  oxygen  in 
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the  catalyst  layer.  In  this  study,  it  is  assumed  that  C qo  is  equal 
the  concentration  in  gas  flow  channel,  fjf1 2  is  exchange  current 
density  at  reference  oxygen  concentration,  and  it  is  calculated 
by  the  following  equation  which  is  obtained  by  Parthasarathy 
etal.  [17]: 


logioO 


3.507 


4001 

T 


(16) 


The  MEA  that  was  used  in  Parthasarathy’s  experiment  was 
different  from  one  of  this  experiment,  so  the  exchange  current 
density  at  reference  oxygen  concentration  of  this  study  is 
calculated  by  the  following  equation: 


:ref+ 

0, 


=  k. 


■ref+ 

‘olO 


(17) 


where  k,0  is  constant,  a.%  of  Eq.  (11)  is  calculated  by  the 
following  equation  [10]: 


a\  =  Ca  +  2.3  x  10“3(T  -  303.15) 


(18) 


Local  current  density  is  obtained  with  the  above  equations. 
But,  in  case  current  density  is  smaller  than  AetQ,  activation 
overvoltage  becomes  negative.  Therefore,  it  has  to  be  taken 
that  this  model  is  not  applicable  on  such  current  density  con¬ 
dition. 

In  this  study,  the  result  and  method  of  Parthasarathy  et  al. 
and  Yoshikawa  et  al.  were  used  in  order  to  calculate  the  activa¬ 
tion  overvoltage.  However,  this  combination  is  inappropriate 
theoretically  because  these  result  and  equation  were  obtained 
in  different  experimental  conditions  from  each  other.  That 
is  why  we  corrected  this  model  by  fitting  several  parame¬ 
ters  with  experimental  data  and  reduced  the  influence  of  the 
gap  between  different  experimental  conditions.  As  the  cur¬ 
rent  density-voltage  curve  shown  in  Figs.  8-11,  it  is  possible 
to  confirm  that  the  influence  of  that  gap  were  low. 


3.  PEFC  thermal  flow  model 

3.1.  Structure  of  analytical  model  and  assumption 

Fig.  4  shows  the  structure  of  PEFC  thermal  flow  model. 
Gas  flow,  concentration  and  temperature  are  calculated  in  the 
gas  flow  channel  of  anode  and  cathode.  Assuming  tempera¬ 
ture  of  membrane,  electrode  and  GDL  are  equal  to  each  other, 
current  density  and  temperature  are  calculated  with  the  com¬ 
bination  of  these  three  components  above.  The  combination 
of  them  is  expressed  to  be  solid  phase.  The  governing  equa¬ 
tions  are  derived  on  the  basis  of  the  following  assumption: 

1 .  The  volume  of  the  condensation  water  is  ignored,  and  the 
water  moves  with  the  gas. 

2.  The  reaction  area  reduction  caused  by  flooding  of  elec¬ 
trode  is  ignored,  and  the  diffusion  prevention  caused  by 
water  condensation  is  also  ignored. 

3.  Fluid  is  incompressible  Newtonian  fluid  and  ideal  gas. 
Flow  condition  is  laminar  flow. 


Cooling 

plate 


Cathode 
gas  phase 
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-*c  •Temperature:  Tb 
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•Velocity: 

vc 
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•  Pressure : 
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Fig.  4.  Schematic  diagram  of  PEFC  reaction  and  flow  analysis  model. 

4.  Heat  transfer  between  separator  and  gas  is  ignored.  But 
heat  transfer  among  gas  phase,  solid  phase  and  back  plate 
is  calculated.  Back  plate  temperature  is  constant. 

5.  Cell  voltage  is  uniform  and  constant. 

6.  Cooling  plates  are  inserted  between  every  three  cells,  and 
the  length  of  heat  conduction  between  cooling  plate  and 
solid  phase  is  lengthened  three  times  as  long  as  the  thick¬ 
ness  of  single  cell. 

7.  The  flux  of  water  through  membrane  is  constant  between 
both  interface  of  anode  side  and  cathode  side. 


3.2.  Derivations  of  governing  equation  of  PEFC  flow 
analysis  model 


The  direction  of  membrane  surface  is  two-dimension 
(x-y),  and  the  two-dimensional  governing  equations  were 
derived.  The  continuity  equation  of  anode  and  cathode  gas  is 
shown  by  the  following  equation: 


(V  •  v)  =  - 


Dp 


(19) 


where  v  is  the  velocity  vector  of  mixture  gas,  /s  the  depth 
of  gas  channel,  p  the  density  of  mixture  gas,  Mj  the  molec¬ 
ular  weight  of  chemical  species  j  and  rj  is  the  reaction  or 
condensation  rate  per  unit  area  of  chemical  species  j. 

The  equation  of  motion  is  shown  by  the  following  expres¬ 
sion: 

=  -S/p  -  [V  •  t]  -  pu(V  ■  v)  (20) 

where  p  is  pressure  and  the  operator  D/Dr  is  substantial  time 
derivative,  r  is  viscous  stress  tensor  that  is  shown  by  the 
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following  expression: 


dvx  2 

Xxx  =  -2^—  +  -m(V  •  v), 
ox  3 

dvy  2 

Xyy  =  -2/iy  +  -M(V  •  1>) 


"txy  —  M 


V  dy  +  3x  ) 


(21) 


where  /i  is  the  viscosity  of  mixture  gas.  And  in  order  to 
include  the  effect  of  viscous  drag  from  both  walls  that  are 
GDL  and  the  bottom  of  channel,  Hele-Shaw  model  shown 
following  was  used: 

_  1 2iivx  drzy  12 fxVy  (22) 

dz  ~  (IS)2  ’  dz  ~  (IS)2 


These  heat  flux  and  heating  value  are  shown  by  the  following 
equation: 


q  3  =  (Eah  —  V)i 

q4  =  h\T&  -  T)  +  hc(Tc  -  Ts) 
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The  equation  of  mass  balance  of  chemical  species  A  is 
shown  by  the  following  equation: 

DCa  i  r  a 

=  T>aV-Ca  -  *  -  CA(V  •  v)  (23) 

where  Ca  is  the  concentration  of  chemical  species  A  and  Da 
is  the  diffusion  coefficient.  The  equations  of  mass  balance 
were  derived  to  eight  kinds  of  C^,  Cj^,  C^20(v),  C^20(1)) 
c02’  Cn2’  Ch20(y)’  CH2od) that  are  hydrogen,  oxygen,  nitro¬ 
gen,  vapor  and  condensed  water  in  anode  and  cathode  chan¬ 
nel. 

The  equation  of  mass  balance  of  condensed  water  is  shown 
by  the  following  equation: 

DCh2o(1)  '-h2o(1)  ^  /V7 

— ^ —  = - ^ - Ch2o(1)(V  ■  v)  (24) 

The  equation  of  energy  of  anode  and  cathode  gas  is  shown 
by  the  following  equation: 

Drg  | 

Pgcpg—  =  fcgv2rg  +  ~(qi  +  q2)  -  pgcgrg(v  •  v)  (25) 

where  Cp  is  the  specific  heat,  T  the  temperature,  k  the  thermal 
conductivity,  the  superscript  g  expresses  anode  gas  or  cathode 
gas,  q\  and  q2  are  heat  flux  from  solid  phase  and  cooling  plate, 
respectively. 


qi  =  h(Ts  -  Tg),  q2  =  h(Tb  -  Tg)  (26) 


where  Eah  is  the  value  of  reduction  change  of  water  enthalpy 
to  voltage,  f'sep  the  heat  conductivity  of  separator,  A  //[[,()  the 
change  of  water  enthalpy  between  vapor  and  liquid  and  bc  is 
the  condensation  rate  constant. 

The  reaction  and  condensation  rates  of  each  ingredient  are 
shown  by  the  following  equation: 
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where  a  is  water  transfer  coefficient. 

Current  density  is  calculated  by  the  above  PEFC  reaction 
model,  which  is  shown  in  Eqs.  ( 1)— ( 18),  and  that  model  is 
shown  by  the  following  equation  with  local  concentration, 
the  temperature  and  saturated  vapor  pressure: 


where  h  is  heat  transfer  coefficient,  7s  the  temperature  of 
solid  phase  and  f°  is  the  temperature  of  cooling  plate. 

The  equation  of  energy  of  solid  phase  is  shown  by  the 
following  equation: 

dTs  1 

PsCp-^-  =  +  —  (<73  +  <74  +  <75  +  qb)  (27) 

where  Is  is  the  thickness  of  solid  phase,  the  superscript  s 
expresses  solid  phase.  q2  is  the  heating  value  owing  to  elec¬ 
trochemical  reaction  per  unit  area  and  time,  q 4  the  heat  flux 
from  gas,  q$  the  heat  flux  from  cooling  plate  and  q^  is  the 
latent  heating  value  of  condensation  per  unit  area  and  time. 


i  =  f(V,  C&, ,  C4o(v),  C^2 ,  rs,  P£2asat)  (30) 

It  is  assumed  that  water  moves  by  electro-osmosis  effect 
and  back  diffusion  effect  in  electrolyte  membrane.  Water 
transfer  coefficient  expresses  the  number  of  net  moving  water 
molecules  from  anode  to  cathode  through  the  membrane, 
when  one  proton  transfers  through  the  membrane.  And  it  is 
calculated  by  the  following  equation  shown  by  Nguyen  and 
White  [5]: 
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where  nd  is  the  electro-osmotic  drag  coefficient,  t^le 

effective  diffusion  of  water  in  membrane,  C^j  the  molar 

concentration  at  cathode  side  interface  of  membrane, 
the  molar  concentration  at  anode  side  interface  of  membrane 
and  tm  is  the  membrane  thickness,  nd  and  DJ")0  are  obtained 
by  the  following  equations: 

nd  =  0.0049  +  2.02§a  -  4.53(fa)2  +  4.09(§a)3  (§a  <  1) 

nd=  1.59  +  0.159(^-1)  (f1  >  1) 

(32) 
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where  D°  is  effective  diffusion  of  water  in  membrane  at  stan¬ 
dard  condition.  C^'q  and  C^'q  are  shown  by  the  following 
equations: 


,A:(m) 

Pi y 

h2o  _ 

tfc(m) 

Ply 

H20  — 
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where  is  the  membrane  dry  density  and  is  the 
membrane  dry  equivalent  weight. 

Eqs.  ( 1)— (34)  are  connected  in  their  respective  group,  such 
as  reaction,  concentration,  temperature  and  flow  and  their  par¬ 
tial  differential  equations  are  discretized  by  finite  differential 
method.  The  boundary  conditions  of  flow  velocity,  tempera¬ 
ture  and  concentration  are  set  as  followings: 

(1)  Gas  inlet  boundary,  these  variables  are  constant. 

(2)  Gas  outlet  boundary :  the  gradients  of  these  variables  are 
constant. 

(3)  Wall  boundary,  the  gradients  of  temperature  and  concen¬ 
tration  are  constant,  and  flow  velocity  is  zero. 

Current  density  and  water  transfer  coefficient  were  calcu¬ 
lated  all  over  the  electrode  area.  In  case  of  calculation  of  these 
variables  at  gas  channel  area,  the  local  concentration  and  tem¬ 
perature  at  that  point  were  used  as  the  needed  variables.  In 
case  of  calculation  of  these  variables  at  no  gas  channel  area 
(at  shoulder  area),  the  local  values  at  adjoining  gas  channel 
were  used.  In  this  case,  the  diffusion  length,  /d,  was  changed  to 

\J  (ld)2  +  (Ax)2  in  Eq.  (4)  that  calculates  the  limiting  current 
density,  and  gas  diffusibility  to  shoulder  area  was  considered. 
Ax  is  the  distance  on  calculation  mesh  to  adjoining  gas  chan¬ 
nel.  And  the  reaction  rate  of  this  shoulder  area  was  put  in  the 
balance  equations  of  adjoining  mesh  point  at  gas  channel, 


Table  3 

Calculation  system  and  method  used  in  this  study 


Analysis  dimension 
Coordinate  system 
Digitization  of  space 
Handing  of  convection  term 
Making  time  dispersed 
Flow  analysis  algorithm 
Iteration  of  pressure 
Mesh 


2-D 

Rectangular  coordinates 
Finite  difference  method 
QUICKER  method 
Adams-Bashforth  scheme 
FS  method 
SOR  method 

Square  staggered  mesh,  480  x  480 


Fig.  5.  Solution  scheme  for  PEFC  analysis. 


and  the  conservation  low  was  satisfied  in  single  cell.  Table  3 
shows  the  calculation  system  and  method  used  in  this  study 
and  Fig.  5  shows  calculation  flow  chart.  Eqs.  (1)— (34)  were 
calculated  until  all  variables  became  constant. 


4.  Water  ejection  experiment 

As  PEFC  is  driven  under  100  °C,  supplied  or  generated 
water  may  be  condensed  and  become  liquid,  and  moreover, 
gas  flow  may  be  inhibited  by  the  condensed  water.  And  it  is 
thought  that  this  problem  results  in  the  cell  instability,  the 
decrease  of  cell  output  and  the  decline  of  durability.  In  this 
study,  water  ejection  condition  was  obtained  by  the  exper¬ 
iment  on  the  assumption  that  condensed  water  block  the 
channel.  Fig.  6  shows  the  experimental  apparatus.  U  shape 
pipe  was  connected  to  the  one  side  of  the  channel  blocked 
by  the  water,  and  length  of  the  blocking  water  was  between 
10  and  20  mm.  Water  drop  was  dripped  to  the  U  shape  pipe, 
and  it  pressurized  the  one  side  of  blocking  water,  and  the 
differential  pressure  that  was  able  to  eject  water  was  mea¬ 
sured.  (This  differential  pressure  is  expressed  as  the  ejection 
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Fig.  6.  Schematic  diagram  of  liquid  water  ejection  experimental  apparatus. 

differential  pressure.)  Three  kinds  of  tube,  which  are  circu¬ 
lar  tube  made  of  halocarbon  resin  (A),  rectangle  tube  made 
of  acrylic  resin  (B)  and  circular  glass  tube  (C),  were  used 
in  this  experiment.  As  a  result,  it  turned  out  that  the  ejec¬ 
tion  differential  pressure  was  not  affected  by  the  blocking 
water  size,  shape  or  material  of  the  tube,  but  by  only  tube 
diameter.  Fig.  7  shows  the  relationship  between  diameter 
and  ejection  differential  pressure.  Diameter  of  rectangular 
tube  was  calculated  as  hydraulic  diameter.  Fig.  7  means  that 
as  the  hydraulic  diameter  decreases,  the  ejection  differential 
pressure  rises.  The  following  ejection  differential  pressure 
equation  was  obtained  in  Fig.  7: 

A P  =  1.04  x  10“2Dh“L42  (35) 

where  A p  is  the  ejection  differential  pressure  and  !\  is  the 
hydraulic  diameter. 


0.1  1  10 
Hydraulic  diameter  [mm] 

Fig.  7.  Relationship  between  the  hydraulic  diameter  and  ejection  differential 
pressure. 


Fig.  8.  Effect  of  O2  concentration  on  cell  voltage/current  density  plots  with 
experimental  and  calculated  results;  cell  temperature  is  333  K,  humidify 
temperature  is  333  K  and  hydrogen  composition  is  100%. 

5.  Results  and  discussions 

5.1.  Experimental  and  calculation  results  ofPEFC 
reaction  model 

In  order  to  construct  the  PEFC  reaction  model,  the  effect  of 
changing  operation  temperature,  humidify  temperature  and 
hydrogen  and  oxygen  concentration  in  supply  gas  on  the  cur¬ 
rent  density  and  voltage  characteristic  of  a  small  PEFC  was 
examined  experimentally.  And  the  parameters  (f  d,fc ,  «a,  «), 
Ca,  Ca  ka,  kj0 )  that  are  applicable  on  all  operation  condition 
were  decided  by  trial  and  error  method  with  experimental 
results  and  Eqs.  ( 1 )— ( 18).  The  decided  values  of  these  param¬ 
eters  are  shown  in  Table  4.  The  diffusion  coefficient  of  hydro- 


Current  density  [A/cm2] 

Fig.  9.  Effect  of  H2  concentration  on  cell  voltage/current  density  plots  with 
experimental  and  calculated  results;  cell  temperature  is  333  K,  humidify 
temperature  is  333  K  and  oxygen  composition  is  21%. 
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Table  4 

Fitting  parameters  decided  by  experiments  with  small  cell 

f“  =  0.020 
fc  =0.043 
aa  =  0.10 
a\  =  0.09 
Ca  =  0.98 
Ca  =  2.0 
^  =  0.15 
ki0  =  1.0 


Fig.  1 0.  Effect  of  cell  temperature  on  cell  voltage/current  density  plots  with 
experimental  and  calculated  results;  humidify  temperature  is  equal  to  cell 
temperature,  hydrogen  composition  is  100%  and  oxygen  composition  is 
21%. 


gen  and  oxygen  were  1.03  x  10-4  and  2.91  x  10-5  m2  s-1, 
respectively,  and  the  standard  electromotive  force  was  1 .23  V. 
Figs.  8  and  9  show  experimental  data  and  calculation  data 
in  case  of  changing  oxygen  and  hydrogen  concentration  in 


Current  density  [A/cm  | 


supply  gas,  respectively.  From  these  graphs,  as  oxygen  and 
hydrogen  concentration  in  supply  gas  were  lower,  the  cell 
voltage  reduced  by  concentration  overvoltage.  Fig.  10  shows 
the  results  in  case  of  changing  cell  and  the  humidifier  tem¬ 
perature.  The  humidify  temperature  was  equal  to  the  cell 
temperature.  In  case  the  cell  temperature  and  the  humidify 
temperature  were  high,  the  cell  voltage  was  high  in  low  cur¬ 
rent  density  because  the  activation  overvoltage  reduced,  but 
the  cell  voltage  was  low  in  high  current  density  because  of  the 
following  cause.  When  the  humidify  temperature  was  high, 
water  vapor  concentration  of  supply  gas  increased,  and  oxy¬ 
gen  concentration  of  supply  gas  decreased  relatively.  And  the 
oxygen  concentration  overvoltage  reduced  the  cell  voltage. 
Fig.  1 1  shows  the  result  of  changing  the  humidifier  temper¬ 
ature.  The  cell  temperature  was  constant.  As  the  humidify 
temperature  were  lower,  the  cell  voltage  reduced  because  of 
low  ionic  conductivity  of  membrane.  In  these  graphs,  the  cal¬ 
culation  results  mostly  agree  with  the  experimental  results. 


Fig.  1 1 .  Effect  of  humidification  temperature  on  cell  voltage/current  density 
plots  with  experimental  and  calculated  results;  cell  temperature  is  363  K, 
hydrogen  composition  is  100%  and  oxygen  composition  is  21%. 


:  The  direction  of  inlet  and  outlet  gas  flow 


Fig.  12.  (A-E)  Shapes  of  five  kinds  of  PEFC  separators. 
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But  hydrogen  concentration  overvoltage  is  discussed  in 
here.  In  this  experiment,  anode  GDL  and  cathode  GDL  were 
equal,  and  it  is  thought  that  anode  effective  porosity,  /a, 
and  cathode  effective  porosity,  /c,  must  be  the  same  essen¬ 
tially.  Moreover,  in  case  that  flooding  by  cathode  generation 
water  strongly  affects  oxygen  diffusion,  it  is  thought  that 
/a  must  be  larger  than  fc.  But  the  repellency  of  the  sur¬ 
face  of  anode  gas  channel  and  GDL  may  have  dropped  in 
this  experiment,  and  so  the  blocking  of  the  hydrogen  diffu¬ 
sion  may  have  occurred.  On  the  other  hand,  it  is  guessed  that 
there  were  the  gap  between  Fick’s  diffusion  model  and  actual 
diffusion  condition,  and  that  hydrogen  diffusion  coefficient 
in  calculation  was  excessive  value.  The  gap  of  calculation 
and  experimental  results  by  changing  hydrogen  concentra¬ 
tion  in  Fig.  9  was  larger  than  that  by  changing  other  condition 
(Figs.  8,  10  and  1 1).  Accordingly,  it  seems  that  the  accuracy 
of  /a  is  lower  than  that  of  other  parameter.  The  next  section 
deals  with  the  calculation  results  in  detail,  and  it  explains  that 
average  hydrogen  concentration  overvoltage  was  about  1% 
of  the  whole  overvoltage,  and  hydrogen  mass  transfer  in  GDL 
did  not  affect  the  whole  cell  system.  So,  it  is  able  to  say  that 
this  problem  about  the  accuracy  of  anode  effective  porosity 
/a  did  not  affect  the  calculation  result  very  much.  However, 


it  is  thought  that  the  anode  concentration  overvoltage  with 
high  accuracy  will  be  needed  in  case  of  other  operating  con¬ 
dition.  In  the  future,  the  accuracy  of  this  model  will  need 
to  be  increased  by  similar  experiments.  Concerning  other 
parameter,  the  PEFC  reaction  model  that  can  estimate  cur¬ 
rent  density  and  voltage  characteristic  on  various  operation 
conditions  was  constructed  with  the  parameters  of  Table  4. 

5.2.  Results  of  PEFC  reaction  and  flow  analysis  in  five 
kinds  of  separator 

Fig.  12  shows  five  kinds  of  separator  shape  that  is  object 
of  this  study.  The  electrode  area  was  a  100  mm  square.  The 
width  of  the  gas  inlet  and  outlet  header  was  5  mm,  and  the 
depth  was  1  mm.  These  separators  were  the  parallel  separator 
A  that  all  channels  were  straight,  serpentine  separator  B  that 
all  channels  were  bunched  and  they  turned  twice,  serpentine 
separator  C  that  all  channels  were  bunched  and  they  turned 
four  times,  semi-serpentine  separator  D  that  was  designed  as 
the  shape  to  combine  the  parallel  separator  with  the  serpen¬ 
tine  separator,  and  the  parallel  separator  E  that  was  partially 
improved  shape  of  separator  A  in  order  to  make  gas  flow 
rate  distribution  uniform.  In  the  concrete,  separator  E  was 


Cathode  separator 


(a) 


Anode  in  Cathode  in 


(b) 


Fig.  13.  Developed  view  of  single  cell  and  gas  flow  pattern  with  separators  A-D.  (a)  Developed  view  of  single  cell  (the  anode  separator  and  the  cathode 
separator  are  the  same);  (b)  gas  flow  pattern  (viewpoint:  from  back  of  the  cathode  separator);  (c)  gas  flow  pattern  (viewpoint:  in  Figs.  15  and  16). 
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Cathode  separator 


(a)  (b) 


Fig.  14.  Developed  view  of  single  cell  and  gas  flow  pattern  with  separator  E.  (a)  Developed  view  of  single  cell  (the  anode  separator  and  the  cathode  separator 
are  different);  (b)  gas  flow  pattern  (viewpoint:  from  back  of  the  cathode  separator);  (c)  gas  flow  pattern  (viewpoint:  in  Figs.  15  and  16). 


designed  to  unify  the  pressure  drop  in  each  channel  by  means 
of  controlling  the  pressure  gradient  in  inlet  and  outlet  header 
section  with  curved  walls.  In  case  of  parallel  separator,  gas 
flow  rate  of  cathode  side  is  larger  than  that  of  anode  side, 
and  it  is  thought  that  the  cathode  gas  flow  rate  in  distant 
channel  from  inlet  is  larger  than  that  in  near  channel  from 
inlet.  Therefore,  the  different  shape  separator  between  anode 
side  and  cathode  side  were  used  on  the  basis  of  such  con¬ 
dition.  Fig.  13  shows  the  developed  view  of  single  cell  and 
gas  flow  pattern  with  separators  A-D.  In  case  of  separators 
A-D,  both  anode  and  cathode  side  separators  were  the  same. 
As  the  same  separators  were  put  together  like  Fig.  13(a),  gas 
flow  patterns  are  shown  in  Fig.  13(b).  And  these  gas  flow 
pattern  from  the  same  point  of  view  in  Figs.  15  and  16  that 
are  current  density  distribution  and  solid  phase  distribution, 
respectively,  are  shown  in  Fig.  13(c).  Accordingly,  in  case  of 
separator  A,  the  gas  flow  pattern  was  parallel  between  anode 
gas  and  cathode  gas  in  reaction  area.  And  Yi  and  Nguyen  [6] 
showed  that  current  density  and  temperature  distributions  in 
case  of  coflow  were  more  uniform  than  those  in  case  of  coun¬ 
terflow.  Fig.  14  shows  the  developed  view  of  single  cell  and 
gas  flow  pattern  with  separator  E.  In  case  of  separator  E,  both 
anode  and  cathode  side  separators  were  different  from  each 


other,  and  the  inlet  and  the  outlet  of  anode  gas  were  differ¬ 
ent  from  those  of  separators  A-D.  As  the  anode  and  cathode 
separators  were  put  together  like  Fig.  14(a),  gas  flow  patterns 
are  shown  in  Fig.  14(b).  As  a  result,  the  gas  flow  pattern  of 
separator  E  was  counter  between  anode  gas  and  cathode  gas 
in  the  reaction  area.  And  this  gas  flow  pattern  from  the  same 
point  of  the  view  in  Figs.  15  and  16  are  shown  in  Fig.  14(c). 
Table  5  shows  the  dimensions  of  MEA,  GDL  and  separator, 
and  Table  6  shows  calculation  condition.  It  was  supposed 
that  physical  properties  were  constant,  because  the  effect  of 
changing  gas  composition  in  cell  on  physical  properties  was 
little. 

Fig.  15  shows  the  calculation  result  of  current  density  dis¬ 
tribution  in  each  separator.  In  this  figure,  the  current  density 
with  separator  A  was  high  at  the  upper  side  of  the  channel 
and  low  at  the  lower  side  of  the  channel.  Because  oxygen 
concentration  were  high  at  the  upper  side  of  the  channel  but 
they  became  low  along  the  channel.  And  current  density  of 
the  center  channel  was  the  lowest,  because  the  gas  flow  rate 
of  the  center  channel  was  the  lowest  in  case  of  separator  A, 
and  the  effect  of  oxygen  concentration  overvoltage  became 
large.  The  current  density  with  separators  B  and  C  was  high 
at  the  center  of  the  cell,  because  this  area  was  middle  of  both 
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Fig.  15.  Calculation  results  of  current  density  distribution. 


Table  5 

Dimensions  of  MEA,  GDL  and  separator  of  single  cell  in  PEFC  reaction 
and  flow  analysis 


MEA 

Thickness  of  membrane  (|xm) 
Size  of  catalyst  layer  (mm2) 
Amount  of  Pt  (g  m-2) 

GDL 

Size  (mm2) 

Thickness  (mm) 

Separator 

Channel  width  (mm) 

Shoulder  width  (mm) 

Channel  depth  (mm) 

Inlet  and  outlet  width  (mm) 
Number  of  channel 

Distance  from  cooling  plate  (mm) 


30 

100  x  100 
3.0 


100  x  100 
300 

1 

1 

1 

5 

A:  50,  B:  16,  C:  10,  D:  16,  E:  50 
7 


anode  and  cathode  gas  channel,  and  anode  relative 
humidity  and  oxygen  concentration  were  not  very  low. 
Therefore,  the  influence  of  the  resistance  and  concentration 
overvoltage  was  little  in  the  area.  As  each  channels  of 
separator  D  meander  separately,  the  area  of  the  high  current 
density  is  close  to  that  of  the  low  current  density,  and  that  is 
distributed.  The  current  density  distribution  with  separator 
E  was  more  uniform  than  that  with  separator  A.  However, 
that  was  not  uniform  perfectly,  and  the  current  density  at  the 
distant  area  from  inlet  was  larger  than  that  at  the  near  area 
from  inlet.  From  the  above-mentioned,  it  followed  that  the 
current  density  distribution  was  different  by  the  difference 
of  the  separator  shape. 

Fig.  16  shows  solid  phase  temperature  distribution  of 
each  separator.  Compared  with  current  density  distribution 
of  Fig.  15,  solid  phase  temperature  was  high  at  the  area 
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Fig.  16.  Temperature  distribution  of  solid  phase. 


of  high  current  density  because  of  the  calorific  value.  In 
this  study,  cooling  plates  were  inserted  between  every  three 
cells,  and  the  length  of  heat  conduction  between  cooling 
plate  and  solid  phase  was  lengthened  three  times  as  long 
as  the  thickness  of  single  cell.  It  was  expected  that  the 
temperature  distribution  of  each  cell  was  more  complex 
by  the  heat  conduction  from  the  neighboring  cell  in  actual 
stack,  but  the  rough  condition  could  be  estimated  by  this 
analysis. 

It  was  thought  that  the  gas  flow  rate  of  cathode  was  not 
more  uniform  than  that  of  anode  because  gas  flow  rate  and 
inertial  force  are  large.  Fig.  17  shows  the  distribution  of 
the  cathode  gas  flow  rate  in  each  separator.  As  the  chan¬ 
nel  number  was  different  in  each  separator,  the  right  edge  of 
a  horizontal  axis  is  expressed  to  the  channel  near  the  inlet. 


the  left  edge  of  a  horizontal  axis  is  expressed  to  the  channel 
near  the  outlet.  The  value  of  vertical  axis  is  each  groove  flow 
rate  divided  by  the  average  flow  rate.  In  case  the  gas  flow 
rate  is  uniform,  this  value  of  vertical  axis  is  1  all  over  the 
channel.  In  the  separator  A,  the  gas  flow  rate  of  the  chan¬ 
nel  that  was  the  nearest  the  outlet  was  3.3  times  as  much 
as  mean  value,  and  the  gas  flow  rate  of  the  central  chan¬ 
nel  was  0.6  times  as  much  as  mean  value.  So,  the  gas  flow 
rate  distribution  of  separator  A  was  the  largest  in  any  sepa¬ 
rators,  and  the  ununiformity  of  the  gas  flow  rate  generated 
the  current  density  distribution  of  separator  A.  Because  of 
improvement  of  header  shape,  the  flow  field  in  separator  E 
is  more  uniform  than  in  separator  A.  But  the  gas  flow  rate 
was  not  uniform  perfectly  even  in  separator  E,  and  the  gas 
flow  rate  of  far  side  from  the  cathode  inlet  was  the  largest. 
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Fig.  17.  Calculation  result  of  gas  flow  rate  distribution. 
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The  gas  flow  rates  of  the  other  separator  were  distributed 
between  0.9  and  1.2  times  as  much  as  mean  value,  and  were 
uniform  comparatively,  because  the  channel  became  longer 
by  meandering  the  channel,  and  pressure  drop  of  the  channel 
became  much  larger  than  that  of  header.  Accordingly,  that 
of  the  separator  C  that  has  the  longest  channel  was  the  most 
uniform. 

The  minimum  pressure  drop  in  the  channel  of  each  separa¬ 
tor  is  shown  in  Fig.  18.  The  horizontal  axis  shows  the  width  of 
the  channel,  and  the  vertical  axis  shows  the  differential  pres¬ 
sure  between  the  inlet  and  the  outlet  of  the  channel.  The  curve 
line  in  Eq.  (35)  obtained  from  the  water  ejection  experiment. 
It  was  found  that  it  was  difficult  to  eject  water  in  separators  A 
and  E  because  the  pressure  drop  was  below  the  ejection  line, 
and  that  it  was  possible  to  eject  it  in  other  separators  because 
the  values  were  above  the  line.  In  addition,  if  the  channel  is 
bent  to  prevent  water  from  blocking,  the  important  thing  is 
not  to  bend  to  excess  but  to  obtain  the  differential  pressure 
exceeding  Eq.  (35). 


Table  6 

Calculation  conditions  of  PEFC  reaction  and  flow  analysis 


Pressure  (MPa) 

0.1 

Operation  temperature  (K) 

333 

Humidifier  temperature  (K) 

333 

Cooling  plate  temperature  (K) 

333 

Operation  voltage  (V) 

0.5 

Inlet  flow  rate  (m3  s  1 ) 

Anode 

1.25  x  10~5 

Cathode 

5.15  x  10“5 

Inlet  gas  composition  (%) 

Anode 

H2:  99.0 

N2:  1.0 

Cathode 

02:  21.0 

N2:  79.0 

Fig.  18.  Relationship  between  separator  shape  and  ejection  of  water. 

5.3.  Comprehensive  evaluation  of  the  performance  of 
five  kinds  of  separator 

Fig.  19  shows  average  current  density  (a),  current  density 
distribution  (b),  solid  phase  temperature  distribution  (c),  gas 
flow  rate  distribution  (d)  and  cathode  pressure  drop  (e)  of 
each  separator.  In  Fig.  19(a),  it  was  found  that  average  current 
densities  were  mutually  equal.  In  Fig.  19(b),  the  values  of  bar 
graph  shows  the  difference  between  the  maximum  value  and 
the  minimum  value  of  current  density,  and  that  of  separator 
A  was  the  largest,  because  the  gas  flow  rate  distribution  was 
more  marked  than  other  separators.  In  Fig.  19(c),  the  values 
show  the  difference  between  the  maximum  value  and  the 
minimum  value  of  the  solid  phase  temperature,  and  that  of 
separator  D  was  the  lowest.  As  the  area  of  high  current  density 
was  close  to  that  of  the  low  current  density  in  separator  D,  the 
temperature  difference  was  reduced  by  the  heat  conduction  in 
the  solid  phase.  On  the  other  hand,  that  of  separators  A  and 
E  was  the  highest,  because  the  heating  value  was  different 
locally  by  ununiform  gas  flow  distribution.  In  Fig.  19(d),  it  is 
a  comparison  of  the  values  in  which  the  maximum  gas  flow 
rate  is  divided  by  minimum  gas  flow  rate  on  cathode,  and 
that  of  separator  A  was  the  largest.  In  Fig.  19(e),  the  cathode 
pressure  drop  of  separator  C  was  the  largest,  and  the  reason 
for  this  was  that  the  channel  of  separator  C  was  longer  than 
that  of  others. 

As  remarkable  current  density  and  temperature  distribu¬ 
tion  have  a  bad  influence  on  the  safety,  durability  and  driving 
efficiency  of  PEFC,  it  is  necessary  to  make  these  distribu¬ 
tions  uniform  as  much  as  possible.  And  in  order  to  make 
the  cell  output  stable,  it  is  important  to  be  able  to  eject  con¬ 
densed  water.  Moreover,  the  pressure  drop  has  to  be  as  low 
as  possible  for  the  purpose  of  reduction  of  gas  supply  power. 
Consequently,  separator  D  that  satisfies  these  necessary  con- 
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(a)  A  B  C  D  E  (b)  A  B  C  D  E 


Fig.  19.  (a-e)  Comprehensive  performance  evaluation  of  five  kinds  of  separators. 


ditions  was  the  best  separator.  However,  the  current  density  6.  Conclusion 
distribution  changes  by  operating  condition,  such  as  cell  tem¬ 
perature,  humidification  and  gas  flow  rate,  so  it  is  necessary  to  In  order  to  make  the  PEFC  reaction  model,  the  effects 

examine  it  in  a  different  condition.  In  addition,  it  is  necessary  of  cell  temperature  and  oxygen  and  hydrogen  concentration 
to  examine  water  ejection  experiment  in  detail  by  changing  on  the  output  characteristic  of  PEFC  were  examined  experi- 

the  material  and  shape  of  channel.  mentally,  and  the  appropriateness  of  calculation  results  was 
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confirmed.  Next,  this  PEFC  reaction  model  was  combined 
with  thermal  flow  analysis  model,  and  PEFC  reaction  and 
flow  analysis  model  of  actual  size  was  made.  Furthermore, 
five  kinds  of  separators  were  evaluated  with  this  model  from 
the  viewpoint  of  gas  flow  condition,  current  density  distribu¬ 
tion,  cell  temperature  distribution,  pressure  drop  and  ejection 
performance  of  liquid  water.  And  the  water  ejection  differen¬ 
tial  pressure  was  investigated  by  experiment.  The  following 
results  were  obtained  by  these  examinations: 

1 .  Average  current  densities  of  each  separator  were  almost 
same. 

2.  The  current  density  distribution  was  caused  by  the  con¬ 
centration  distribution  resulted  from  ununiformity  of  gas 
flow  rate. 

3.  In  parallel  channel  separator,  as  gas  flow  rate  distribution 
was  remarkable,  the  distributions  of  concentration,  cur¬ 
rent  density  and  solid  phase  temperature  became  large, 
and  it  is  difficult  to  eject  liquid  water.  But  those  distri¬ 
butions  were  decreased  by  improvement  of  the  shape  of 
header. 

4.  In  serpentine  channel  separator,  gas  flow  rate,  current 
density  and  temperature  were  more  uniform  than  that  of 
parallel  separator.  But,  even  if  the  channel  is  lengthened, 
it  does  not  influence  these  distributions  so  much,  and  the 
pressure  drop  increases. 

5.  In  semi-serpentine  channel  separator,  the  parts  of  high  and 
low  current  density  were  dispersed  and  close  to  each  other. 
And  the  solid  phase  temperature  was  unified. 

The  PEFC  reaction  model  was  developed  for  the  pur¬ 
pose  of  simulating  actual  cell  condition  by  fitting  sev¬ 
eral  parameters  with  experimental  results.  But  this  model 
has  many  assumptions.  In  our  future  study,  it  is  expected 
that  this  model  will  be  improved  to  apply  in  the  various 
conditions. 
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